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It was recently suggested that molecular junctions would be excellent elements for efficient and 
high-power thermoelectric energy conversion devices. However, experimental measurements of ther¬ 
moelectric conversion in molecular junctions have indicated rather poor efficiency, raising the ques¬ 
tion of whether it is indeed possible to design a setup for molecular junctions that will exhibit 
enhanced thermoelectric performance. Here we suggest that hybrid single-molecule nanoparticle 
junctions can serve as efficient thermoelectric converters. The introduction of a semiconducting 
nanoparticle introduces new tuning capabilities, which are absent in conventional metal-molecule- 
metal junctions. Using a generic model for the molecule and nanoparticle with realistic parameters, 
we demonstrate that the thermopower can be of the order of hundreds of microvolts per degree 
Kelvin, and that the thermoelectric figure of merit can reach values close to one, an improvement 
of four orders of magnitude improvement over existing measurements. This favorable performance 
persists over a wide range of experimentally relevant parameters, and is robust against disorder 
(in the form of surface-attached molecules) and against electron decoherence at the nanoparticle 
molecule interface. 


INTRODUCTION 

In recent years, the study of single-molecule junctions 
- the ultimate limit of electronic nano-teclmology has 
progressed well beyond their role as functional elements 
in electronic devices [1], and today additional function¬ 
alities, from opto-electronics and spintronics, through 
phononics, to thermoelectricity [2-4] are under investiga¬ 
tion. The potential applicability of single-molecule junc¬ 
tions in thermoelectricity - the conversion of heat into 
electric power - is of particular interest, since thermo¬ 
electricity may turn out to be an important element in 
addressing global energy issues [5, 6]. The clear advan¬ 
tages of thermoelectric energy conversion, such as the 
’’green” nature of the energy conversion process, its ap¬ 
plicability for waste heat harvesting, and the easy device 
maintenance due to the absence of moving parts, would 
lead us to think that thermoelectric devices should al¬ 
ready be a substantial part of the energy market, yet this 
is not the case. The reason lies in the simple fact that 
current thermoelectric devices are not efficient enough, 
making competition with traditional (large scale) energy 
conversion systems virtually impossible. 

Due to their versatility, low dimensionality, and low 
thermal conductivity, molecular junctions offer a possi¬ 
ble route for enhancing thermoelectric performance [7]. 
Indeed, various theoretical studies have suggested that 
molecular junctions could be a key component of efficient 
and high-power thermoelectric devices [8-16]. Experi¬ 
mental demonstrations [17 26], however, have indicated 
just the opposite, as can be deduced by examining the 
two central parameters of thermoelectric conversion, the 
thermopower (or the Seebeck coefficient) and the ther¬ 
moelectric figure of merit (FOM). 


The thermopower S measures the voltage generated 
per unit temperature difference (in the linear response) 
[4, 19, 27, 28]. Values of S ~ 10 2 — 10 3 /jV /K are typical 
for standard semi-conductor based thermoelectrics, yet 
molecular junctions exhibit small values of S , typically 
S ~ 5 — 50 fiV/ K. The FOM, namely ZT, defined as 
ZT = ( ffj T , where G is the conductance of the junction, 
k = K e + Kphn is the total thermal conductance, which 
includes both electronic (e) and phononic (phn) contri¬ 
butions, and T is the temperature. ZT is directly related 
to the device efficiency [29], and ZT —» oo corresponds 
to the Carnot efficiency (thus, theoretically, there is no 
upper bound on ZT). It is commonly held that, from the 
efficiency perspective, a ZT of ~ 4 is required for thermo¬ 
electric conversion to be competitive [30]. However, typ¬ 
ical ZT values obtained from measurements in molecular 
junctions are ZT ~ 10 -3 — 10~ 5 . Thus, there seems to 
be a discrepancy between theoretical and computational 
studies of thermoelectric conversion in molecular junc¬ 
tions, which in many cases [8-14] predict S ~ 10 2 — 10 3 
/jV/K and FOM of ZT 1, and the experimental eval¬ 
uation of the two measures. 

The origin of this discrepancy seems to stem from two 
factors. First, one has to be careful to take the phonon 
contribution to the thermal conductance into account; 
a typical and realistic value for the phonon thermal con¬ 
ductance in molecular junctions is K p h n = 10—100 pW/K 
[13, 19, 31 -33]. Second, many calculations show en¬ 
hanced thermoelectric performance based on tuning of 
the molecular orbitals or, equivalently, the Fermi level 
of the electrodes. However, in reality this tuning is very 
difficult to achieve [24, 34-36], and the only ’’tuning pa¬ 
rameter” available for molecular junctions is the choice 
of the molecular moiety. It is thus a central challenge to 
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design molecular junctions that are both tunable in some 
way and show favorable thermoelectric performance with 
realistic parameters. 

In this paper, we present a new setup for molecular 
junctions that has the potential to achieve this goal. 
Our setup is based on a hybrid single-molecule semi¬ 
conducting nanoparticle (SC-NP) structure in which a 
molecule is connected on one side to a metallic electrode 
(as in conventional molecular junctions), and on the other 
side connected to the second electrode through a SC-NP 
placed on the electrode, as is schematically depicted in 
Fig. 1(a). In this setup, the molecular junction can be 
tuned via tuning the electronic properties of the nanopar¬ 
ticle, namely, by choosing a suitable material and by con¬ 
trolling the size and shape of the nanoparticle (for reviews 
see, e.g., [37-39]). Fabrication of hybrid NP-single moel- 
cule junctions have already been demonstrated with Au 
NPs (see, e.g., [40-42]), making our suggested system 
experimentally feasible. 

Using a generic model for the SC-NP, we show that this 
junction can reach values of S ~ 100 — 400 /zV/K and 
ZT ~ 1 for a broad range of realistic parameters. The 
origin of the enhanced thermoelectric performance can be 
traced to the interplay between the local transport prop¬ 
erties of the molecule and the gapped density of states 
(DOS) of the nanoparticle [13]. Thermoelectricity typi¬ 
cally requires a large particle-hole asymmetry (reflected 
in the fact that, at low temperatures, the thermopower is 
proportional to the derivative of the transmission func¬ 
tion [4, 19, 28]). This asymmetry is enhanced in this 
junction due to the presence of the semiconducting gap 
in the SC-NP [13], an effect that is rectified due to the 
finite size of the nanoparticle. We show that the opti¬ 
mal parameters for thermoelectric conversion depend on 
the geometry of the nanoparticle and the contact geome¬ 
try between the nanoparticle and the molecule. Further, 
we demonstrate that the favorable thermoelectric per¬ 
formance is robust against disorder (in the form of sur¬ 
face dangling molecules) and dephasing, and finally, we 
discuss the temperature dependence of the FOM, which 
exhibits a maximum at T ~ 450 K. 


MODEL AND CALCULATION 

The transport and thermoelectric properties of the 
hybrid molecule-nanoparticle junction are calculated by 
using the non-equilibrium Green’s function approach, 
which has become the standard tool for such calculations 
[27, 28]. The junction (graphically depicted in Fig. 1(a)) 
is described using the Hamiltonian: 

P = Pb+Pnp+Pm+Pt+Pb-np+Pnp-m+Pm-t , 

( 1 ) 

that includes the bottom electrode (B), the nanoparticle 
(NP), the molecule (M), the top electrode (T), and the 


coupling between the bottom electrode and nanoparticle 
(B-NP), between the NP and the molecule (NP-M) and 
between the molecule and the top electrode (M-T). Since 
electron spin does not play a role in the mechanisms we 
describe here for enhancement of thermoelectricity, we 
treat spinless electrons. 

To describe the SC-NP, we use a generic tight-binding 
model for a semiconductor [43-48] of the form 

Pnp = + y( _1 ) Pr )4cr - t ^2 clc r > , (2) 

r (rr'> 

where the summation is taken over the atom positions r 
(assumed to form a cubic lattice); cj(c r ) creates (annihi¬ 
lates) an electron at position r; t is the nearest-neighbor 
hopping matrix element; e c is the position of the band- 
gap center; and A is the semiconductor band-gap. The 
function P r = x + y + z gives a modulation of the on-site 
energy e c ± 4 between neighboring atoms, generating a 
gapped band at the thermodynamic limit [44]. We point 
that although ab initio methods for calculating proper¬ 
ties of NPs are emerging [49, 50], these are still not de¬ 
veloped for transport calculations. Furthermore, since 
we are aiming at presenting general properties of hybrid 
junctions, our tight-binding calculation is generic and not 
limited to a specific system. 

The molecule is described as a single orbital, which can 
correspond to either the HOMO or LUMO, depending 
on the position of the orbital energy [28]. The molecular 
Hamiltonian is simply 

Pm = co d^d , (3) 

where d f ( d ) creates (annihilates) an electron at the 
molecular orbital. In addition to the calculations de¬ 
scribed below, we performed calculations that also in¬ 
clude a Coulomb interaction term (for the molecule) 
and spin-full fermions (by using the equations-of-motion 
method [51, 52]) but found no quantitative change in 
the main results, and therefore we keep the description 
here as simple as possible and treat electrons as non¬ 
interacting particles. The coupling between the molecule 
and the nanoparticle is described by the Hamiltonian: 

Pnp-m = —t\c\ M d + h.c. , (4) 

where ym is the position of the atom in the nanoparticle 
that is in contact with the molecule, and t\ is the hopping 
matrix element related to the overlap integral between 
the molecular orbital and the atomic level at r m- 

The metallic top and bottom electrodes are assumed 
to be non-interacting metals [53] with the Hamilto¬ 
nian P x = J2k e k c l,x c k,x, X = T,B [27, 28, 53]. 
The coupling between the molecule and the top elec¬ 
trode (representing, e.g., the tip of a scanning tunnel¬ 
ing microscope) is represented by the Hamiltonian term 
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Hm-t = Efc t V k,Tc{ T d + h.c.. Similarly, the con¬ 
tact between the bottom electrode and the nanoparti¬ 
cle is described by Hb-mp = E r -es ^ k,Bc' kB c r + h.c.. 
To proceed with the transport calculations, the metallic 
electrodes are treated within a wide-band approximation 
[54] (valid since our model is for non-interacting elec¬ 
trodes) , where the top electrode is defined by a (retarded 
and advanced) self-energy term £^ a = |, 

in which | M) is the single-particle molecular orbital, 
and Ft is the electrode-induced level broadening. Simi¬ 
larly, the bottom electrode is defined via the self-energy 
term £^° = =Ryf- Eree |r)(r|, where T# is the broad¬ 
ening due to the bottom electrode and the summation 
is taken over all the atom positions in the nanoparti¬ 
cle that are in contact with the bottom electrode. We 
point that the model described above (and specifically 
the wide band approximation) implies electron-phonon- 
induced thermalization in the electrodes, but no electron- 
phonon interaction in the junction at this stage. All the 
relevant energies- the Fermi level of the electrodes, NP 
valence and conduction bands, and molecular orbitals - 
are schematically also shown in Fig. 1(a). 

Once the Hamiltonian and the self-energies are de¬ 
fined, the calculation proceeds via the non-equilibrium 
Green’s function approach, which is reduced to the Lan- 
dauer formalism for non-interacting electrons. [27, 55, 
56]. The Green’s functions are determined via G r,a = 
(. E-n + £ r ’°) _1 , where £ r ’ a = £j a + £^ a . The trans¬ 
mission function is given by T(E) = Tr (£^G r £^G°), 
and the transport coefficients, namely, the conductance 
G, the thermopower S, and the thermal conductance 
ft, are determined within the Landauer formalism as 


G = e 2 Loi S = Li/(eTLo, k = (^L 2 — / T , where 

T is the temperature (room temperature, unless other¬ 
wise stated) and L n = — ^ f dET(E)(E — /z)”^g are the 
Landauer integrals, with h being Planck’s constant, /z 
the chemical potential of the electrodes, and f{E) the 
Fermi-Dirac distributions. The thermoelectric FOM ZT 
is given by ZT = ^. 


RESULTS 

We start by describing the thermopower S and FOM 
ZT for a single molecule placed on a square pyramid¬ 
shaped nanoparticle [57, 58] (the bottom electrode is in 
contact with the [111] plane), as shown in Fig. 1(a). Since 
we present here a generic model for nanoparticles, we are 
not aiming at obtaining quantitative results describing 
a specific system. We are, nonetheless, aware that it is 
essential to take numerical parameters that are realis¬ 
tic and readily describe experimental systems. We thus 
choose the semiconducting band center at e c = —4.8 eV 
and A = 0.8 eV, corresponding to PbSe nanoparticles, 
and /z = —5.1 eV as the electrode chemical potential 


(corresponding to Au electrodes). Other numerical pa¬ 
rameters were t = 1.6 eV , T# = 0.05 eV, Tt = 0.01 
eV (describing weakly coupled molecules, see, e.g.,[59]). 
The pyramid basis contains 10 x 10 atoms, and our re¬ 
sults only weakly (and quantitatively) depend on the size 
of the nanoparticle. 

In Fig. 1(b), we show the thermopower A as a function 
of molecular orbital energy level eg and the molecule- 
nanoparticle coupling t\. In experiments, eo can be tuned 
by choice of molecule and by choice of the nanoparticle 
composition and size. The coupling t\ can be addition¬ 
ally tuned by stretching or squeezing the molecular junc¬ 
tion with the top electrode [60-62], As may be seen, S 
can reach values as high as ±200 /zV/K and can change 
sign according to the position of the molecular level with 
respect to the semiconducting band edge. 

The appearance of a thermopower maximum upon a 
change in eo is not surprising, since one would expect that 
tuning eo would lead to a near resonance in transmission 
(which implies a maximum thermopower). However, the 
appearance of a maximum upon a change of t\ is sur¬ 
prising; in a typical single-molecule junction, stretching 
the junction will only lead to a change in the molecule- 
electrode coupling, and will nor result in a thermopower 
maximum. Here, since the molecule and the NP hy¬ 
bridize, changing t± (experimentally - by means of pulling 
the junction) is similar to changing the molecular orbital, 
i.e. pulling the junction plays the role of gating. Since 
gating a molecular junction is a challenging task [24, 34- 
36], this result puts an additional advantage on hybrid 
junctions. 

To calculate ZT, we add to the electronic thermal con¬ 
ductance a phononic term K p h n = 50 pW/K, a realistic 
value for molecular junctions [19, 31-33]. In Fig. 2, we 
plot the central result of this paper, ZT of the hybrid 
single-molecule-NP junction, as a function of eo and t±. 
The FOM reaches ZT > 0.8, i.e., more than three orders 
of magnitude larger than the values measured in regular 
molecular junctions. The realistic parameters considered 
here and the wide range of parameters in which ZT is 
large implies that this regime should be accessible for 
experiments. 

In the inset of Fig. 2, we show ZT as a function of eo for 
a constant t\ = 0.1 eV for three different configurations. 
The first (solid blue line) is the same as the setup in the 
main figure. In the second setup (dotted orange line), the 
molecule is in contact with a pyramid nanoparticle, but 
the positions of the ±A terms in the nanoparticle Hamil¬ 
tonian of Eq. 2 are switched, modeling a change in the 
atom species at the apex of the nanoparticle (for instance 
either Pb or Se in PbSe nanoparticles). The third setup 
(dashed green line) describes a molecule in contact with a 
cube-shaped (as opposed to a pyramid) nanoparticle. As 
may be seen, the contact configuration and the shape of 
the nanoparticle can have a strong effect on ZT ; although 
values of ZT > 0.8 can be achieved in these configura- 
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FIG. 1. (a) Schematic representation of the hybrid metal- 
molecule-nanoparticle-metal junctions (the choice of molecule 
is arbitrary), and the corresponding energy landscape, (b) 
Thermopower S as a function of the position of the molecular 
orbital eo and molecule-nanoparticle coupling fi for a hybrid 
molecular junction with a square-pyramid-shape (see text for 
numerical parameters). 


tions, the optimal parameters vary between the different 
setups. 



FIG. 2. ZT as a function of the position of the molecular or¬ 
bital eo and the molecule-nanoparticle coupling fi for a hybrid 
molecular junction with a square-pyramid shape (see text for 
numerical parameters). Inset: ZT as a function of eo for a 
constant ti =0.1 eV, for the hybrid junction with a pyramid 
nanoparticle (blue solid line), a pyramid nanoparticle with 
switched ±A positions (dotted orange line, see text), and a 
cube-shaped nanoparticle (dashed green line). 


In experiments, a reasonable situation would be for ad¬ 


ditional molecules to attach to the surface of the nanopar¬ 
ticle. To model this scenario, we add to the Hamiltonian 
additional molecules, with the same orbital level eo and 
coupling ti but with a random coupling point to the sur¬ 
face of the nanoparticle (top inset in Fig. 3). The trans¬ 
port properties are then averaged over 10 4 realizations of 
random positions. In Fig. 3, ZT is plotted as a function 
of surface coverage (in percents) of attached molecules. 
Surprisingly, we found a slight increase of ZT for a small 
number (~ 8%) of attached molecules, followed by a de¬ 
crease in ZT when the number of attached molecules was 
increased. To elucidate the origin of this result, the lower 
inset shows the average conductance (blue circles) and 
thermopower (orange triangles) as a function of surface 
coverage (in percents). We found that the average con¬ 
ductance actually increased with the number of surface 
molecules, but the thermopower decreased, thus eventu¬ 
ally leading to a decrease in ZT. 

The origin of this effect is the fact that the surface cov¬ 
erage induces two competing processes. On one hand, the 
addition of molecules bound to the surface adds conduc¬ 
tion channels (i.e. local resonances in the transmission 
function) and therefore increases the conductance. One 
the other hand, the presence of disorder tends to flat¬ 
ten the resonances on average, and as a result the ther¬ 
mopower, which is proportional to the derivative of the 
transmission function [4], and hence becomes smaller as 
the transmission resonance becomes wider, is reduced. 
This competition is reflected in the opposite trends of G 
and S in the inset of Fig. 3. Since ZT is the product of 
the conductance which increases and the thermopower 
which decreases, it exhibits non-monotonic behavior. 



surface coverage [%] 


FIG. 3. ZT as a function of surface coverage (in percents) of 
molecules attached to the nanoparticle surface. Upper inset: 
schematic representation of the surface-attached molecules. 
Lower inset: average conductance (blue circles) and ther¬ 
mopower (orange triangles) as a function of number of surface 
molecules. 

At room temperature it is possible that the electron 
motion across the junction will not be coherent and that 
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the electrons will dephase as they cross the molecule-NP 
interface. This may occur due to the interaction of elec¬ 
trons with the (soft) phonons of the NP, or due to vibra¬ 
tions in the position of the molecule with respect to the 
electrodes (but not due to interaction of electrons with 
the molecular vibrations, which are typically high-energy 
modes; such interactions were discussed in, e.g. [63]). If 
this is the case, the formulation presented above is in¬ 
valid, as it describes coherent transport. To account for 
decoherence, we note that if the electron loses its phase 
on the molecule-NP interface, then the NP can, in fact, 
be considered as a semiconducting electrode. This re¬ 
sults in an effective SC-molecule-metal junction, which 
can again be described using the Landauer formula, with 
the appropriate choice of self-energies. 

To model the SC electrode, we recall that the imagi¬ 
nary part of the self-energy describes the electrode DOS. 
Similar to the wide-bancl approximation for the metal¬ 
lic electrode, we model the semiconducting electrode as 
constant outside the gap and zero in the gap. The imag¬ 
inary part of the semiconducting self-energy can thus be 
written with the use of a lieaviside step function O as 
IrnSsc = yf- (0(|I? — e c | — y)), and the real part of the 
self-energy determined via the Kramers-Kronig relation 
[27, 64]. In the numerical calculations given below, the 
step-function discontinuity is broadened by 0.01 eV, a 
value that can arise in realistic systems from lattice im¬ 
purities and dislocations or thermal fluctuations. 

In Fig. 4 we plot the (a) conductance, (b) ther¬ 
mopower, (c) thermal conductance and (d) ZT, for the 
hybrid SC-NP-molecule-metal junction (solid lines) as 
a function of e 0 , with Tb = 0.05 eV, Tt = 0.01 eV, 
e c = —4.8 eV, and A = 0.8 eV. Looking at ZT, we find 
that the dephasing process in the NP does not substan¬ 
tially reduce the maximal ZT from the coherent case (the 
maximal value of Fig. 2). For comparison, the dashed 
lines of Fig. 4(a-d) show the same for A = 0, i.e. a ’’stan¬ 
dard” metal-molecule-metal (M-M-M) junction. The two 
most striking features of the comparison between the hy¬ 
brid junction and the M-M-M junction are (i) the ther¬ 
mopower is substantially larger for the hybrid junction 


As was noted earlier, one of the advantages of the hy¬ 
brid NP-molecule junctions is the ability to tune not 
the molecular orbitals, but the SC-NP band structure. 
With this in mind, in Fig. 5 we plot ZT for a constant 
value of the molecular orbital eo = —5.2 eV, as a func¬ 
tion of the band-gap A (the band center is at e c = —4.8 
eV and A = 0 corresponds to a M-M-M junction, see 
Fig. 1(a)), exploring the two cases of weakly coupled 
molecule (solid blue line) and strongly coupled molecule 
(dashed red line), as described above. For the weakly 
coupled molecule we see an increase of a factor 2 from 


and over a wider range of eo, and (ii) ZT of the hybrid 
junction is roughly twice as large as that of the M-M-M 
junction. 

The inset of Fig. 4(d) shows the inverse Lorenz number 
y as a function of eo for the hybrid junction (solid line) 
and M-M-M junction (dashed line). While a violation 
of the Wiedeman-Franz (WF) law ^ ~ 1 is observed 
for both types of junctions, the hybrid junction shows a 
much larger violation, which in fact determines the posi¬ 
tion of the maximal ZT, while for the M-M-M junction, 
the position of the maximal ZT is determined by the 
maximum of S. 

Up till now, we considered a molecule which is weakly 
coupled to the electrodes. However, depending on the 
chemical moiety, the coupling between the molecule and 
the electrodes can be much larger (see, e.g. [9, 10, 59]). 
It is therefore of interest to see whether the increase in 
ZT in hybrid junctions compared to the M-M-M junc¬ 
tions is maintained also for strongly coupled molecules. 
In Fig. 4(e-h), the same as in Fig. 4(a-d) is plotted for a 
strongly coupled molecule, with Tb =0.5 eV, Tt = 0.1 
eV. Again, the two most striking differences between the 
hybrid junction and the M-M-M junction are the ther¬ 
mopower and ZT, which are even more profound for the 
strongly coupled junction. Due to the large value of the 
coupling, the transmission resonance of the M-M-M junc¬ 
tion is very broad, and therefore the thermopower is very 
small. That, in addition with the fact that the WF law 
is obeyed (shown in the inset of Fig. 4(h), dashed line) 
results in a small ZT ~ 10~ 2 . In contrast, the hybrid 
junction (solid lines) shows only slight reduction of the 
thermopower (Fig. 4(f)), because although the molecular 
level is broadened due to the metallic lead, the band-edge 
of the SC electrode still introduces a sharp feature to the 
transmission function. Along with a violation of the WF 
law (inset of Fig. 4(h)), this leads to a relatively large 
ZT ~ 1. Surprisingly, ZT for the strongly-coupled hy¬ 
brid junction displays larger ZT than the weakly-coupled 
junction, with large values of ZT for eo well inside the 
SC band gap (as a result of the WF law violation inside 
the gap). 


the M-M-M junction to the optimal band-gap. For the 
strongly coupled junction we find that ZT can rise as 
high as ~ 0.6, with a four orders of magnitude increase in 
ZT compared to the M-M-M junctions. The inset shows 
S' as a function of A, and both the weakly-coupled and 
strongly-coupled molecules exhibit orders-of-magnitude 
increase in S compared to the M-M-M junctions. 

Finally, in Fig. 6, the temperature dependence of ZT is 
examined for the weakly-coupled hybrid molecular junc¬ 
tion, evaluated for eo = —5.2 eV (corresponding to the 
maximum in ZT from Fig. 4(d)) (the rest of the pa- 
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FIG. 4. (a) conductance, (b) thermopower, (c) thermal conductance and (d) ZT as a function of molecular orbital energy eo, 
for a semiconductor-molecule-metal junction (solid lines) and a metal-molecule-metal junction (dashed lines), calculated for a 
weakly coupled molecule (Tg = 0.05 eV, Ft = 0.01 eV). Inset of (d): Inverse Lorenz number as a function of eo. (e-h) same as 
(a-d), for a strongly coupled molecule (Tb = 0.05 eV, IV = 0.01 eV). 



A [eV] 


FIG. 5. ZT as a function of the SC NP band-gap A at con¬ 
stant molecular level eo = —5.2 eV, for weakly coupled (solid 
blue line) and strongly coupled (dashed red line) molecule. 
Inset: S as a function of A. 


rameters are the same as in Fig. 4(a-d)). ZT is found 
to increases with temperature, exhibiting a maximum 
of ZT ~ 2 at T ~ 450 K, followed by a moderate de¬ 
crease. This finding again implies that relatively large 
values of ZT persist in a broad range of parameters. In 
fact, in the inset of Fig. 5, we plot ZT calculated with an 
over-evaluated value of the phonon thermal conductance 
Kphn = 150 pW/K, and we find that ZT reaches values 
of ZT ~ 0.3, which is still several orders of magnitude 
larger than the observed values. 



HKj 


FIG. 6. Temperature dependence of ZT for eo = —5.16 eV. 
Inset: same, but with K v hn = 150 pW/K. 


SUMMARY AND CONCLUSIONS 

In summary, we presented calculations of thermopower 
and thermoelectric FOM for a hybrid metal-single 
molecule- semiconducting nanoparticle-metal molecular 
junction. The presence of the SC-NP and the position of 
the molecular orbital close to the semiconducting band 
edge enhanced the particle-hole asymmetry required for 
efficient thermoelectric conversion. The resulting values 
of thermopower and ZT were much larger than those 
measured in experiments, and reach values as high as 
S ~ 500 /iV/K and ZT ~ 1 (an improvement of four 
orders of magnitude improvement over measured molec¬ 
ular junctions). We showed that this enhanced thermo¬ 
electric performance persists over a wide range of pa¬ 
rameters and is robust against disorder in the form of 
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surface-attached molecules. We showed that decoherence 
at the nrolecule-nanoparticle boundary is not detrimental 
to the thermoelectric performance, and that large values 
of ZT persist up to high temperatures. Comparing hy¬ 
brid SC-NP-molecule junctions to the more ’’standard” 
metal-molecule-metal junction, we found that for weakly- 
coupled molecules there is a factor ~ 2 increase in ZT and 
a factor ~ 4 in the thermopower. For strongly coupled 
molecules, the advantage of hybrid NP-molecule junc¬ 
tions is even more profound, with 2 — 4 orders of magni¬ 
tude increase in ZT and thermopower n hybrid junctions. 

The model we presented here is a generic model, not 
aimed at any specific system. Nevertheless, our numeri¬ 
cal parameters were taken from experimentally observed 
value, including the phononic contribution to the thermal 
conductance. This, along with the fact that enhanced 
thermoelectric performance was found for a wide range 
of molecular parameters and was robust against disor¬ 
der, decoherence and high temperatures, strongly sug¬ 
gest that high values of thermopower can be reached in 
future experiments on nrolecule-nanoparticle junctions, 
which are promising candidates for nano-scale thermo¬ 
electric conversion. 
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